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The University of Houston-Clear Lake established the Research Institute for 
Computing and Information systems in 1986 to encourage NASA Johnson Space 
Center and local industry to actively support research in the computing and 
information sciences. As part of this endeavor, UH-Clear Lake proposed a 
partnership with JSC to jointly define and manage an integrated program of research 
in advanced data processing technology needed for JSC’s main missions, including 
administrative, engineering and science responsibilities. JSC agreed and entered into 
a three-year cooperative agreement with UH-Clear Lake beginning in May, 1 986, to 
jointly plan and execute such research through RICIS. Additionally, under 
Cooperative Agreement NCC 9-16, computing and educational facilities are shared 
by the two institutions to conduct the research - - 

The mission of RICIS is to conduct, coordinate and disseminate research oh 
computing and information systems among researchers, sponsors and users from 
UH-Clear Lake, NASA/JSC, and other research organizations. Within UH-Clear 
Lake, the mission is being implemented through interdisciplinary involvement of 
faculty and students from each of the four schools: Business, Education, Human 
Sciences and Humanities, and Natural and Applied Sciences. 

Other research organizations are involved via the “gateway” concept. UH-Clear 
Lake establishes relationships with other universities and research organ izationsT 
having common research interests, to provide additional sources of expertise to 
conduct needed research. 

A major role of RICIS is to find the best match of sponsors, researchers and 
research objectives to advance knowledge in the computing and information^ 
sciences. Working jointly with NASA/JSC, RICIS advises on research needs, 
recommends principals for conducting the research, provides technical and 
administrative support to coordinate the research, and integrates technical results 
into the cooperative goals of UH-Clear Lake and NAS A / JSC. 
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Abstract 

Necessary and sufficient conditions for positive realness in 
terms of state space matrices are presented under the assumption of 
complete controllability and complete observability of square 
systems with independent inputs. As an alternative to the positive 
real lemma and to the s-domain inequalities, these conditions 
provide a recursive algorithm for testing positive realness which 
result in a set of simple algebraic conditions. By relating the 
positive real property to the associated variational problem, the 
paper outlines a unified derivation of necessary and sufficient 
conditions for optimality of both singular and nonsingular problems. 

1. Introduction 

Positive real systems play a major role in control theory, 
especially in adaptive control, and in stability analysis. The 
impressive development of adaptive control and self-turning 
regulation over the last two decades {1,2] is hinged on satisfaction 
of some positive realness conditions. Alternatively, considerable 
initial knowledge about the controlled plant must be given. The 
prior knowledge is used to implement reference models, identifiers, 
or observer-based controllers of about the same order as the plant 
Since the prior assumptions about the controlled plant may never be 
entirely satisfied, the stability properties of the related adaptive 
schemes are debatable. Therefore, a direct adaptive control 
procedure which docs not use identifier or observer-based 
controllers in the feedback loop is preferred The implementation of 
such an algorithm requires positive real controlled plants or 
alternatively, a synthesis of a positive real plant on the basis of the 
actual plant 


1 . This work was partially sponsored by Eglin AFB under Contract F08635-87* 
K04 17 and by NASA Johnson Space Center through the R1CIS Program of the 
University of Housoo at Gear Lake. 

2. Department of Aerospace Engineering and Engineering Mechanics, University 
of Texas. Austin. TX 78712. 

3. Mechanical, Aerospace and Nuclear Engineering Department, University of 
California. Los Angeles. CA 90024-1597 

4. RAFAEL P. O B 2250, Haifa. Isrtal 


The existing tools for analysis and synthesis of positive real 
systems are based in the s-domain on complex variable inequalities 
which are inconvenient or in the state space requiring the positive 
real lemma equations. These tools are computationally complex and 
there i$ a need for an easily used complementary tool. In Sections 2 
and 3, necessary and sufficient conditions for positive real systems 
with independent inputs are developed using optimal control theory 
for the associated partially singular problem. It is shown that in the 
totally singular case, these conditions are consistent with the 
generalized Legendre -Gebsch condition {3,4]. The new conditions 
are associated with the state space matrices of a minimal realization 
of a square system. The resulting test for positive realness reduces 
to recursively testing certain square matrices for positive definiteness 
and the solution to an algebraic Riccati equation. As an immediate 
result of the new necessary and sufficient conditions, we also show 
that the zeros of a positive real system lie in the closed left half 
complex plane. Some examples are given in Section 4 to illustrate 
the theory. Concluding marks are given in Section 5. 

The derivation of the above results is relaied to dissipative 
systems. Basic definitions and physical characteristics are presented 
below. 

1.1 Dissipative System 

Consider the system input-output description H: U -» Y 
where U * Lj (R+) and Y * LJJ (R+). The notation \\ (R„) is used 
to denote the space of square imegrable functions f: R„ — ► R' 
where R*. * (tp,~). The supply rate associated with this system is 
defined as a function w: R / x R* -* R where 

w(u.y) * y'Qy + 2y'Su + u*Ru (1.1) 

and Q € R**" , Se R"**. Re R w are constant matrices, with 
Q and R symmetric. 

Definition 1.1 [5]: A dynamical system H is dissipative 
with respect to the supply rate w(u,y) if and only if 

" 

I w[u(t), y(t)] dt > 0 (1.2) 

io 

for all ti > to and all u e L^, whenever the initial state satisfies 
x(to) = 0. 





Remark 1. 1: Passivity corresponds to dissipativeness 
where Q = R = 0. 1 ■= m. S * j I m and is mxm identity matrix. 

Remark 1.2: Positive realness corresponds to passivity 
where the dynamical system is linear and time invariant 

Remark U: The concept of a supply rate is related in the 
general case to the "stored energy" for the system. As an example* 
suppose that the system under consideration is an electrical network* 
whose elements are constants, and y(t) the vector of corresponding 
port voltage. Then the system is dissipative with respect to the 
supply rate w(u,y) * u*y provided that all the resistances, 
inductances and capacitance are non-negative. 

1.2 Energy, Power and Information Relationships in 

Dissipative Systems 

The class of dissipative systems which has a finite 
dimensional Internal state is completely described in terms of energy 
storage and power dissipation. Considering this class, the various 
facets of the standard state space model can be associated with the 
concepts of energy, power and information. 

Assume that the system under consideration Is described by 
a linear, time-invariant system 

**Ax + Bu (13) 

y*Cx + Du (1.4) 

where x e R* , u e R ; , y c R" and A, B, C and D are constant 
matrices with appropriate dimension. Then, following [6], the 
system matrices can be regarded as repres enting: 

1. an energy-transformation and dissipation map, associated 
with the matrix A. 

2. a power injection map, associated with the matrices B 

and D. 

3. an information -ex traction map, associated with the 
matrix C 

Figure 1 describes the energy-power-information maps 
associated with the system matrices. 



Fig. 1 Encrgy-Power-Information maps associated 
with the System Matrices 

The matrix B represents the input coupling between the 
information represented by the applied input signals and the power 
available for injection into the system states. The matrix C 
represents the output coupling between the energy in the system 
states and the information in the available output signals. The matrix 
D represents the output coupling between the information 


represented by the applied input signals and the injected power into 
the available output signals. 

1 . 3 Review of the Positive Real Property 

The positive real property is related directly to the transfer 
function matrix description of the system. The positive real lemma, 
presented in Section 2, connects the positive realness to the 
parameters of a system realization with complete controllability and 
complete observability. 

The Positive Real Property J7J: Let G(s) be an m x m 
matrix of functions of a complex variable s. then G(s) is termed 
positive real if the following conditions are satisfied: 

(i) All the elements of G(s) are analytic in Re[s] > 0. 

(ii) G(s) is real for real positive s. 

(iii) G*(s) + G(s) £ 0 for Re [s] > 0 , 
where (*)* denotes complex conjugate transpose. 

Remark 1.4: If G(s) is a real rational matrix of functions 
of s, then necessary and sufficient conditions for the positive real 
property to hold are given by the following theorem. ^ 

Theorem LI {7J: Let G($) be a real rational matrix of 
functions of s. Then, G(s) is positive real if and only if: 

(i) No element of G(s) has a pole in Re{s] > 0. 

(ii) G*(j<a) + G(ja>) £ 0 for all real co, with jcu not a pole of 
any element of G($). 

(iii) If j<Uo ts a pole of any element of G(s), it is at most a 
simple pole, and the residue matrix. 


lira 

S->j<A0 

(s-jcoo)G(s) 

if jcoo is finite. 

lira 

G(s)/s 

if joio is infinite. 


is nonnegative definite Hermitian, 

Following Definition 1.1, if the system is positive real, the 
angle between the output vector y(t) and the input vector u(t) is 
bounded below by - 90 deg. and above by ♦ 90 deg. 

2, Relations Between Optimal Control 
and Positive JR eal ness 

2,1 The Related Variational Problem 
Consider the cost functional 

Vlxo,k>,u(0] = j ! w[u(t). y(t)] dl (2. 1 ) 

to 

where the supply rate 

w (u,y) » y'u = u'D'u + x'Cu (2.2) 

is associated with system (1.3) and (1.4), where the dimensions of 
u and y are m. The problem is to find necessary and sufficient 
conditions for optimality of u* (•) e U to minimize V(xo,io.u( )|, 






denoted V*[xo,to], subject to the dynamic equation of (1.3) where 
X(to) = xo is prescribed. 

Remark 2.1: Since only the symmetric part of D 

contributes 10 w(u,y). then 

w(u.y) - j ( u'Ru + 2*'Cu ) (2.3) 


where 


R = D + D' (2.4) 

Remark 2.2: If R £ 0* and rank (R) = r < m, there exists 
an orthogonal transformation f = [ T, , Tj] such that 



R[Ti.r 2 ] 



(2.5) 


where R r is positive. For instance* T, and Tj may consist of 
normalized eigenvectors of R associated with nonzero and zero 
eigenvalues, respectively [8]. There is a natural partitioning of the 
control vector associated with this transformation* a r-dimensional 
nonsingular control and an (m-r><ii mention al singular control 


2.2 Positive Real Lemma Equations 

Necessary and sufficient condition for V*[xo,to) to be 
bounded below over a finite time interval [to. ti] are presented in 
Theorem IL3.3 of [9]. The required positive real conditions are 
obtained via the extension of the optimality condition to the time- 
invariant, infinite -time case [10]. 

Under the complete controllability and complete 
observability assumption of system (1.3), necessary and sufficient 
conditions for the nonnegativity of V[ 0, to. u(»)l are that there exist 
x < 0, L, and W such that 


xA + A’n 
B*x + C 


jcB + C' 
R 


-[^[UWliO. 


( 2 - 6 ) 


where W and L are matrices with proper dimension. 

By identifying P * - x, the positive real Lemma is stated. 
The Positive Real Lemma (7): Let G(s) be an mxra 
matrix of real rational functions of a complex variable s, with 
G(«o) < ••. Let (A* B* C* D) be a minimal realization of G(s). 
Then, G(s) is positive real if and only if there exist real matrices P, 
L, and W with P positive definite and symmetric, such than 


PA + AP « - L/L 

(2.7) 

BT> > C - WL 

(2.8) 

WW -D + D' 

(2.9) 


Remark 2 2: The generalized Legendre -Oebsch condition, 
which is a necessary condition for V*[xo,tol > - 00 in the totally 
singular case, given in [3] for a linear rime-invariant system can be 


written as 

^(Hu) = CB-(CB)‘ = 0 (2.10) 

|j(Hu)-CAB + (CAB)'S0 (2.11) 

where H is the variational Hamiltonian and X € R" is the associated 


Lagrange multiplier 

H = u'Cx + X'( Ax + Bu ), X' = -H„. 

By letting R « 0. the necessary conditions (2.10) and (2.1 1) are 
also obtained from the positive real lemma. 


3. Positive Real Conditions in Terms 

of State-Space Matrices 

* 

Necessary and sufficient conditions for the nonnegativity of 
V[ 0, to, u( )] are given by the existence of x < 0, L, and W which 
satisfy (2.6). Let G(s) be an mxm matrix of degree n. Consider a 
minimal realization (A, B, C, D) representing the finite-dimensional 
linear time-invariant dynamic equations given by (1.3) and ( 1 .4). In 
terms of state space matrices A, B, C, and D, (2.6) gives necessary 
and sufficient conditions for a positive real system. In this section, 
new necessary and sufficient conditions are developed. 


3.1 Standard Formulation of the Partially Singular 
Problem 

Assume that G(s) is a square matrix of proper rational 
function with independent columns. For any realization, the 
matrices C and B are full rank. Without loss of generality, wc 
consider a minimal realization { A, B, C, D) of the form that A is an 

nxn matrix which is partitioned as 

f A|| Aj2 "I 

L A^j A 22 J ’ 

where An is a kxk matrix, and A 22 is an (n-k)x(n-k) matrix, 


where k is the dimension of the singular control, 

fCrl 


B = [ B r , B* ], C 


-El- 


and 


D + D' 



where R r is a rxr nonsingular matrix corresponding to the 
nonsingular control, B r is an nxr matrix, B t is an nxk matrix 
related to the singular control, Q is a rxn matrix, and C* is a kxn 


matrix, where r = m-k is the dimension of the nonsingular control. 
If n > k, then C, has the following form 
C # -[C,i, 0], 

where C»i is a nonsingular matrix. Correspondingly, B, is written 



We define this as a standard realization. 


Notice that the realization can be obtained by choosing 


suitable bases for the state space and the input/oufput space. For 
example, suppose {A, B, C, G(~)) is a minimal realization of 
G(s). Let the column vectors of f, where f is described in 


Remark 2.2, be a basis of the input/output space, then the 
following transformation y = Tq, u = Tv is defined. Furthermore, 

let qi. q 2 qa-k. qn-k*i* ... , qn be a basis of the state space, 

where q n -k*i span the null space of r i 'c, *nd qi • q 2 

q n .k are arbitrary vectors such that Q = [qj, q 2 , — . qn k* q n -k*t. 


qn ] is nonsingular. This defines a transformation x = Q%. The 
resulting dynamic equation can be written as 


^ = A^ + Bv (3.1) 

q = + Dv, (3.2) 


where A = Q ! AQ, B = Q ‘BP. C = TCQ. and D = PG( ~>r. The 
transfer function matrix of this system is PG(s)f. the positive 


i 



realness of G(s) is equivalent to the positive realness of rG(s)r. 
The application of (2.6) and development of the new necessary and 
sufficient conditions for the partially singular problem will be 
discussed under assumption of a standard realization as discussed. 

3.2 Derivation of New Necessary and Sufficient 
Conditions 

Necessary and sufficient condition for nonnegative of 
V(0,t<>.u(*)] as given by condition (2.6) can be restated in the 


I ' ' 

’ *11 

X 

”1 

I 0 

0 I 

. *t 2* 

*1 J 

. -(ri)-'ki2’ I 



0] 

L 0 

*1 J 


L 

0 Xt J 


following equivalent forms: There exist a x < 0 and a matrix V xB + C * 

X u B*i + * 120*2 


r c.r 


' 0 ‘ 

such that 

. Xn'Bjj + x 5 B * 2 


L 0 . 


A 


■ V'V. 


(3.3) 


'V/V f . 


(3.6) 


it A + A 'n xB + C* 

B'lt + C R 
Furthermore, R being positive semi-definite is a necessary 
condition for satisfying (3.3). If R > 0, then (3.3) can be reduced 
to a condition based upon a Reccari equation. That is, these exists a 
negative definite solution it to the algebraic Riccad equation 
ii( A - BR* ] C ) + ( A - BR^C Vic - x BR-iB*x - CR^C = 0 . 

(3.4) 

If R is singular, (3.3) can be written as 

xA + A'x xB r + C r ' zB,+ C,' 

B r 'x + C f R* 0 * V'V 

L B,’x ♦ C* 0 0 

or, equivalently, there exist a it < 0 and a matrix V r such that 

xB, + C/ -0 (3.5) 

and 

[ itA + A'x JtB r + Cj 
Bf x 4 - Cf R* 

If the dimension of the state is leu than or equal to the 
dimension of the singular control, Le„ nSk, x can be determined 
from equation (3.5). If and only if a x < 0 is solvable from (3.5) 
and the same x satisfies (3.6), the system is positive reaL If n > 
k, the fact x < 0 and equation (3.5) imply that 

C,B, - (C*B,)' - - B*’x B s > 0. (3.7) 

Since C, ■ [ Cn, 0 ], and Qj is nonsingular. Equation (3.7) 
also implies that B,j is nonsingular. Furthermore, (3.5) provides a 
linear constraint on x which is in Lemma 3.1 below. 

Lemma 3.1: x<0, x B f ♦ C*’ - 0 If and only if C*B, 

> Oand 

% m [ -(B.i') 1 C i |+(B f i')-«B u , jc I B.j<B i ,)- 1 •(B„) iB,j , )t, 1 

L (B,j)-> *, J 

(3.8) 

To prove 

(3.9) 

(3.10) 


for some Xj < 0 . 

Proof: Denote x as 


*n 

* 12 * 


*12 

*1 


sufficiency, we assume that X| < 0, C*B, > 0, and 

*u * • (B.O ’C.i + (B*i')' l B* 2 ’X]B* 2 (B $ i )- 1 
*12 = * (B,|’)* I B I 2 , X|. 

I -X 12 (X|)** 

0 I 


Define F = 


FxF = 


Since C,B, * C*iB,i > 0, 

^ (B,jT ! C,i » - (B.iT'C.iB, (B # j )- 1 < 0 . 
Therefore FxF < 0, and it also implies that x < 0. Furthermore, 
by using Xji and Xjj defined in (3.9) and (3.10), we get 


(3.11) 


Next, we prove the necessity. If x < 0, then Xi < 0. From xB, + 
C,' aOwe get 

*nB|i + x u B l2 + C t \' = 0 (3.12) 

* 12 B*i + X!B l2 =*0 (3.13) 

By solving (3.12) and (3.13), the expressions of x n and x !2 are 

obtained which are the same as shown in equations ( 3 . 9 ) and 
(3.10). 

Q.E.D. 


MU.R*): 


Let the matrix shown in (3.6) be denoted as M( x, R r ) 
xA + A'x xB f ♦ C r ‘ 

B r 'x + C r Rf 
For any nonsingular matrix T, (3.6) is equivalent to TM( x, R r )T 

■ Vy Vt, where Vy is a matrix with proper dimension. By defining 
'0 B,i 0] 

T r « I B*j 0 

.001 

and using x defined in (3.8) as a function of xj, then T r is 
nonsingular, and 

f M U M t2 Mj 2 
Ml 2 M 22 M 23 

M 13* M 2 3 ' M 33 J 


T/M(x, R f )T, 


where 


Mu - £ 0 I 0 ] M(x,R) 


-0 

I 

.oj 


1 *i( A 22 - B, 2 (B*i)‘ 1 Aij) ♦ (A 22 - B,2(B s i )* 1 A i2)'xi 


M 12 -[0 I 0 ] M(x,R) 


B»i 
B»2 
0 J 


r *1 ( A 2 iB,i+ A 22 B i2 - 


Bi 2 (B f i)' l A|iB,i-B* 2 <B l |)* 1 Ai 2 B f2 ) * C*iAt 2 


M22 * [ B*i\ B*2*» 0 ] M(x,R) 
- - (C*AB, + B,*A'Ci*) 

M23 - [ Bif, B, 2 \ 0 ] M(x,R) 
« . C,B r ♦ B 


, then F is nonsingular and 


M 33 s [ 0. 0, 1 ] M(x.R) 


01 

0 


B*i 
B*2 
L 0 


ro 

0 

1 J 


Rr- 





— 

By defining 



Ai * A 22 * B^CBgi^An 

(3.14) 


Bi = [ A2iB tl + A22B»2- B t 2(B I i)* I AiiB l f 


— 

B^gO^AnBgj, Bfz] 

(3.15) 


Ci =[-Q,Ai2. 0] 

(3.16) 


R| 


1. 


(C,AB,+ B.'A'C,') -C.B.+B.’C,' 

- B,'C,'+C r B, R, 

a condition which is equivalent to (3.6) can be stated as the follows: 
There exist a jcj <0 and a matrix V! such that 

JljA 1 + Aj'Kj ttiBr+Cj* 

B i'Jt| + Ci Rt 


VfV,. 


(3.17) 


(3.18) 


According to the positive real lemma. Equation (3.18) implies that 
R 

( Ai.BuCi, is positive real. 


A i = A 22 * B, 2 (B,i) i A 1 2 is equal to the zero polynomial of the 
system up to a nonzero scalar factor. 

w(s) 

Proof: Let det G(s) = det ( C(Is - A) >B) = — . *here 

A(s) 


A(s) = det ( Is - A ) 

and y(s) is the zero polynomial of the system. Since state feedbacks 
do not change the numerator of the transfer function matrix, for any 
matrix K, 

det (G k (s)) = det ( C(Is - A - BK )>B) = 

A k (s) 


where A k (s) = det ( Is - A - BK ) 

Let K = ( 0, (B* l )' 1 A I2 l, then 

An 0 
A 21 Aj . 


A + BK 


J An 

. A 71 


A22 - 12 J 




3.3 Necessary and Sufficient Conditions for Positive 

Realness 

The results in Section 3.2 are summarized in the next 
theorem as an alternative necessary and sufficient condition for 
testing positive realness of a square system. 

Theorem 3.1: The necessary and sufficient condition for 
{A, B, C, D) to be positive real is that 

(i) R £ 0; 

(ii) If R > 0* there exists a positive definite solution P to the 
following algebraic Riccati equation 

P( A - BR^C ) + ( A- BR* l C )‘P + P BR ^BT + CR *C - 0; 

(iii) If rank R * r < m, andnSm-r, there exists 

P = CjBgXBjBs’)" 1 > 0 satisfying PB f * <Y and 
f - P A - A'P - PB r + C/ 1 

L - B,’P + C, R, ] i0; 

(iv) If rank R = r < m, and n > m-r, then C*B, ■ (C,B,)* > 

R 

0 and { A|, Bi, C|, ) is positive real, where Aj, Bi, Ci, and 

Ri arc defined in equations (3.14) to (3.17). 

Condition (ii) is obtained by identifying P with - ft in 
equation (3.4). Condition (iii) is the interpretation of (3.5) and (3.6) 
for the case n £ (m - r). If P ® - ft > 0 exists, then PB S » C*\ 
PBgB,' * Cg'B,, and P * CgBgXBgB^)* 1 > 0 . Condition (iv) 
corresponds to the situation we discussed through (3.7) to (3. 13) 

Remark 3.1: Alternative transformation approaches to the 
singular problem using the Kelley transformation for the linear 
quadratic problem are given in [9] for the matrix case. The approach 
here is different via the structure of ft given by Lemma 3.1. 

Remark 3.2: If {A, B. C. D) is a minimal realization, 
then it is required for a positive real system that there exists a 
positive definite matrix P such that 

PA + AP £ 0 

Therefore, it is required that Re [A] £ 0 and the Jordan form of 
A has no blocks of size greater than Ixl with pure imaginary 
diagonal elements. 

Remark 3.3: If G(s) is strictly proper, the minimal 
realization is totally singular, then the characteristic polynomial of 


A k (s) * det ( Is - A - BK ) * det 


Is-Au 0 
-A 2 i Is-A, . 


= det ( Is- An) det (Is-Aj) 


det (G k (s)) = det ( C(b - A - BK ) 'B) « det (C,,(Is-A n ) •' B u ) 
det (C, t ) det (B,i) 
det (Is-Au) 


Therefore, 

y(s) « A k (s) det (G k (s)) * det (Qj) det (B fl ) det (Is-Aj) 

Q.E.D. 


Remark 3.4: From (3.18) and Remark 3.3, we conclude 
that there are n - m finite zeros for a positive real system and all the 
zeros lie in the closed left half complex plane. In other words, the 
system is minimum phase. 


4. Examples 


Theorem 3.1 introduces a recursive procedure for testing 
positive real systems, requests only for testing a series of matrices 
C it Bu > 0, for i = 0, 1, 2, ... , 1, and the solution to a algebraic 
Riccati equation Pi > 0, where i is the index associated with the 
new system obtained from the i-th iteration, and i = 0 corresponds 
loBg, C,, andP. The testing stops when R| becomes non singular, 
or the dimension of the state is less or equal to the dimension of the 
singular control. 

The following examples illustrate the application of Theorem 
3.1. 

Example 4.1: Given G(s) = s ( /+ f)fsT 3) 1 an 


observable realization of G(s) is 



1 

l 

at 

0 
1 


T 

A = 

-3 0 1 

. B = 

4 


0 0 0 


.4. 


Fust iteration: 

R =0 

CB = (CB) f = 1 > 0 


l -ST) 





Ri = 0. 


5. Summary and Conclusions 


■f - j 1 1. 

L- 4 oj 


Second iteration: 

Rj *0 

C|B,«(C,B ,)*--! <0 

Therefore, the system is noc positive real. 


Example 4.2: Given G(s) = ^(f-f 4) • “ 

observable realization of G(s) is 



" -6 1 0 ’ 


T 

A = 

-8 0 1 

, B« 

2 


.000- 


.1. 


C » [ 1,0,0], D * 0. 


First iteranon: 

R = 0 

CB = (CB)’ = 1 > 0 


-2 1 

r It 


At- 

b HJ 

. C, =.[-1,0], R, = 8. 

-1 0 

L4J 


Second iteration: 


-15 

8 1 

‘ -15 

-4 “ 

P| 

+ 


-4 

0 J 

8 

0 a 


Pi + 


Pi + 


1 4 

4 16 

which has a positive definite solution 

Pi 


1 0 
0 0 


0.0394 -0.0225 


l -0.225 0.1557 J 
Therefore, the system is positive real. 

S 2 + I* 


>0 


Example 43: 

realization of G(s) is 


G(s) 


s (S 2 + p2) * 


A minimal 



0 1 

0 “ 


-r 

A 3 

-p 2 0 

1 

. B- 

0 


0 

0 

0 „ 


.z 2 . 


A, = 


’ Bl= [ 2 0 p2 ]- 


C - 1 1. 0. 0 ], D = 0. 


First iteration: 

R *0 

CB « (CB)' * 1 > 0 

0 1 

-z 2 0 

Ri e 0. 

Second iteration: 

Ri =0. 

CiBi * ( CiBi)’ = p 2 -z 2 > 0 if and only if p^>z 2 
A 2 ■ 0, B 2 - • z 2 (z 2 -p 2 ), C 2 = 1, R 2 « 0. 

Third iteration: 

R 2 = 0. 

Pj= C 2 B2 , (B 2 B 2 , ) > = -f ( p2-z2 - > 0 if p J > z 2 
- P2A2 - A2 P2 = 0. 

Therefore, the system is positive real if and only if p 2 > z 2 


This paper reviews positive real system as a subclass of 
dissipative systems and states the positive real lemma equations. By 
using the variational problem associated with the partially singular 
problem, necessary and sufficient conditions for a system to be 
positive real are derived. These conditions are panicul arly 
transparent by using Lemma 3.1 which provides a uniquely 
structure for the matrix x. These positive realness conditions arc 
expressed in terms of the state space matrix inequalities and 
algebraic Riccati equations and do not deal with inequalities in the s 
domain or with solutions of the positive real lemma equations. 
These tests are direct, and a system either satisfies these conditions 
or noL There is no requirement to search over all matrices to 
determine if a condition can be satisfied as in the positive real 
lemma. Examples are given which demonstrate the power of this 
approach. 
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